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ABSTRACT 

We present the first results of a program to characterize the disk and envelope structure of typical 
Class protostars in nearby low-mass star forming regions. We use Spitzer IRS mid-infrared spectra, 
high resolution CARMA 230 GHz continuum imaging, and 2-D radiative transfer models to constrain 
the envelope structure, as well as the size and mass of the circum-protostellar disk in Serpens FIRS 1. 
The primary envelope parameters (centrifugal radius, outer radius, outflow opening angle, and incli- 
nation) are well constrained by the spectral energy distribution (SED), including Spitzer IRAC and 
MIPS photometry, IRS spectra, and 1.1 mm Bolocam photometry. These together with the excel- 
lent Mw-coverage (4.5 — 500 kA) of multiple antenna configurations with CARMA allow for a robust 
separation of the envelope and a resolved disk. The SED of Serpens FIRS 1 is best fit by an enve- 
lope with the density profile of a rotating, collapsing spheroid with an inner (centrifugal) radius of 
approximately 600 AU, and the millimeter data by a large resolved disk with Mdisk ~ 1.0 M© and 
Rdisk ~ 300 AU. These results suggest that large, massive disks can be present early in the main 
accretion phase. Results for the larger, unbiased sample of Class sources in the Perseus, Serpens, 
and Ophiuchus molecular clouds are needed to determine if relatively massive disks are typical in the 
Class stage. 

Subject headings: stars: formation — ISM: individual (Serpens) - submillimeter - infrared: ISM 



1. INTRODUCTION 

Protostars build up their mass by accreting material 
from a dense protostellar envelope, presumably via a ro- 
tationally supported circum-protostellar accretion disk. 
Disk formation is a natural result of collapse in a ro- 
tating core, but it is not known how soon after pro- 
tostellar formation the disk appears, or how massive 
it is at early times. Theory suggests that centrifu- 
gally supported disks should start out small (radius less 
than 10 AU), and thus l ow in ass, and grow with time 
(|Terebev. Shu, fc Casseiil 11984( 1. Unstable or magneti- 
cally supported disks, however, could be much larger 
(radii up to 1000 AU in the magnetically supported case; 
iGalli fc Shulll993l ). and thus more massive at early times. 

The remnants of these protostellar accretion disks are 
easily observed in more evolved phases (e.g. TTauri 
stars), but given that the majority of mass is accreted 
during earlier embedded phases, understanding disks at 
early times is critical. Directly observing disks during 
this main accretion phase is quite difficult, however, as 
they are hidden within dense, extincting protostellar en- 
velopes. The structure of the envelope at small radii is 
another important characteristic of main accretion phase 
protostars that is similarly difficult to directly observe. 
Disk growth or the presence of a binary companion may 
clear out the inner region of the envelope early on, as 
infe rred for the binary Cl ass source IRAS 16293-2422 
by ()j0rgensen et al.ll2005l ). 

There has been a recent push toward detecting disks 



in more embedded objects, with many now known 
and roughly characterized in Class I protostar s (e.g. 
Loo nev ct al. 2000; J0rgensen et al. 2005b; Eis ner et al.l 
20051: lAndrews fc Williamsl [2007D and a few detected 
in the earlier Class stage (e.g. iChandler et all 119951: 
Harvev et all 120031 iBrown et all 120001 : iLoonev et all 
200011 . ' Most previous detailed studies have been lim- 
ited to the most well-known or brightest Class sources, 
however, due to instrumental limitations and a lack of 
large unbiased target samples. The ong oing Submillime- 
ter Array s urvey of low-mass protostars (jJorgensen et al] 
HMSIlooi) is a notable exception. 

With recent large surveys of nearby molecular clouds 
at mid-infrared and (sub)millimeter wavelengths it is 
now possible to define complete samples of Class 
protostar s based on lumino s ity or envelope mass 
limit s (e.g. iHatcheh et al.l l2007t iJorgensen et "all l200a 
iDunham et al.l l200a lEnoch et all 120091 : lEvans et all 

We have recently begun a campaign to characterize 
disk properties in a large, uniform sample of Class 
protostars in nearby low mass star-forming regions (M. 

^ We use definitions of Class 0, Class I, an d Class II ([Andre i 
1 19941) based on t he bolometric temperature (M yers fc Laddlll993l : 
IChen et al.|[T995l) : T,,,,; < 70 K (Class 0); 70 k< rj,^, < 650 K 
(Class I); 650 K< Tf,oi < 2800 K (Class II). We further as- 
sume that classes co rrespond to an evolutionary sequence (e.g. 
IRobitaille et al.l 120061 ): in Class the protostar has accreted less 
than half its final mass (Af« < M^nv), in Class I M« > Menv, and 
in Class II the envelope has dispersed, leaving only a circum-stellar 
disk. 
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Fig. 1. — Spitzer 24 /im image of the immediate environment of Serpens FIRS 1, in the Serpens main core. Bolocam 1.1 mm continuum 
contours are overlaid. Submillimeter source designations for the brightest Casali ct al. (1993) sources are indicated. The nearest embedded 
protostar to FIRS 1 is approx imately 45" or 11000 AU away (Ser-emb 12: lEnoch et al.ii2009l '). and the nearest YSO is 25" or 6000 AU away 
(c2d 142: IHarvev et a.U\^mh . 



Enoch etal. 2009, in preparation). Our study is based on 
the complete (to envelope masses > 0.1 Mq) sample of 
39 Class protostars in the Serpens, Perseus, and Ophi - 
uchus molecular clouds, identified bv lEnoch et al.l ()2009D 
by comparing large-scale Spitzer IRAC and MIPS and 
Bolocam 1.1 mm continuum surveys of the three clouds. 

Combining Spitzer IRS mid-infrared (MIR) spectra 
and high resolution CARMA 230 GHz continuum imag- 
ing with radiative transfer modeling of this sample will 
help to address several fundamental questions about the 
structure and evolution of the youngest protostars: 1) 
How soon after the initial collapse of the parent core does 
a circum-protostellar disk form? 2) What fraction of the 
total circum-protostellar mass resides in the disk, and 
does this fraction vary with time? 3) Are large "holes" 
in the inner envelope, such as that found for IRAS 16293 
bv'Jorgen sen et all (|2005f ). common at early times? 

MIR spectra and millimeter maps provide complemen- 
tary approaches to these questions. The amount of flux 
escaping at A < 50 /^m from deeply embedded sources is 
very sensitive to the opacity clos e to the protostar, an d 
thus the envelope structure (e.g. |j0rgensen et al."2005D. 
While the MIR flux is insensitive to disk properties, high 
resolution millimeter continuum mapping can directly 
detect emission from dust grains in the disk. Millimeter 
observations with excellent M?;-coverage, combined with 
radiative transfer models, can separate the disk from the 
envelope and constrain the disk mass and size. 

Our ultimate goal is to characterize the disk mass, size, 
and inner envelope structure of typical low-mass Class 
protostars, and to quantify any trends with evolutionary 



indicators. In this initial paper we present results for 
Serpens FIRS 1, a well known Class source, which will 
serve as a test case for the full program. 

2. SERPENS FIRS 1 

FIRS 1 is located at 18''29"49.6^ + 01°15'21.9" 
(J2000)^ in the main core (Cluster A) of the Serpens 
Molecular Clotid. We adopt a distance of d = 260 ± 10 pc 
(jStraizvs et al.l 119960 . and any quoted literature values 
are scaled to this distance. It is a well know Class 
protostar (e.g. Hurt fc Barsonv'l l996r) fi rst noted in the 
far- infrared bv IHarvevT W ilking. fc Jov l ([l984), and also 
known by its sub-m illimeter designation Serpens SMM 1 
(|Casah et al.|[T99l . There is a narrow A — 3.6 cm bipo- 
lar radio jet at the position of FIRS 1 ([Rodriguez et al.l 
119891: iCuriel. Moran. fc Cantol 11993'). indicating a pow- 
erful outflow that is also clearly seen in molecular lines 
CDavi s et al.|[T999l: ICuriel et al.| ]l996) . 

Figure [T] gives an overview of the FIRS 1 environment 
with Spitzer 24 fim and Bolocam 1.1 mm continuum 
maps of the Serpens main core. The nearest known YSO 
is appr oximately 25'^ away , or 6000 AU in projected dis- 
tance ("H arvey et al.l [20071 ). and the nearest embedded 
protostar known to have an envelope is 45" « 11000 AU 
away (Eno ch et aI]|2009D . 

FIR S 1 is referred to as Ser-emb 6 in lEnoch et al.l 
([2009[ ). and is a ssociated with th e 1.1 mm Bolocam 
core Ser-Bolo 23 ([Enoch et al.ll2007l) . Based on 2MASS, 
Spitzer, and Bolocam data the bolometric luminosity is 

Spitzer position from (Harvey et al.l II2007I ) 
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Fig. 2.— Spitzer IRS spectrum of Serpens FIRS 1, using the Low Res 7.4 - 14.5 ^im (SLl), Hi Res 9.9 - 19.6 ^im (SH), and Hi Res 
18.7 — 37.2 /im (LH) modules. The low signal-to-noise SH data at A < 13.5 fira is not plotted. Binned data (AA ~ 1 fim) are over-plotted 
as diamonds; error bars represent the error in the mean for each bin. Binned fluxes are used in the model fitting and given in Table [T] 

11.0 Lq,'^ the bolometric temperature is 56 K, confirm- 
ing the Class designation, and the total envelope mass 
is 8.OM0 (Enoch ct al. 2009). 

Previous high resolution millimeter observations have 
placed limits on the mass of a compact disk in FIRS 1. 
iHogerheijde. van Dishoeck. fc Salverdal (|1999f l used ob- 
servation from the Owens Valley Radio Observatory 
(OVRO) millimeter interferometer to estimate a total 
mass (disk p lus envelope) within 100 AU of 0.7 Mq. 
iBrown et al.l ([200Q) place a lower limit on the disk mass 
of ~ 0.1 Mq with sub- millimeter observations {v ~ 
350 GHz) from the James Clerk Maxwell Telescope- 
Cahech Submilhmeter Observatory (JCMT-CSO) single 
baseline interferometer. 

3. OBSERVATIONS 

3.1. Spitzer IRS spectrum 

Mid-infrared spectra were obtained with the Infrared 
Spectrograph (IRS ) on the Sp itzer Space Telescope 
(jWerner et al.ll2004t IHouck et al.|[200l during 2007 Oc- 
tober 5 with the Low Res 7.4-14.5 fim (SLl; R ~ 100), 
Hi Res 9.9 - 19.6 /zm (SH; R - 600), and Hi Res 
18.7 - 37.2 ^m (LH; R ~ 600) modules. Integration 
times were 117 sec in SLl, 189 sec in SH, and 59 sec 

^ Note tha t low er resolution HIRES IRAS fluxes 
l|Hurt fc Barsonvl 119961) and ISO-LWS spectra C Larsson et al.l 
120001) vield higher bolometric luminosities of 62 Lq and 95 Lq, 
respectively, but these may be confused with nearby embedded 
protostars. 



in LH. Off-source or background spectra with the same 
integration times were also obtained for SH and LH. 

Spectra were extracted from the SSC pipeline ver- 
sion S16.1.0 BCD images using the reduction pipeline 
(jLahuis et all 120061 ) developed for the "From Molecular 
Cores to Planet-forming Disks" Spitzer Leg acy Program 
("Cores to Disks" or c2d: lEvans et al.ll2003l i. We use the 
optimal PSF extraction method of the pipeline, which 
is based on fitting the analytical cross dispersion point 
spread function, plus extended emission, interpolating 
over bad pixels. The 1-D spectra are flux calibrated us- 
ing a spectral response function derived from a suite of 
calibrator stars, corrected for instrumental fringe resid- 
uals, and an empirical order matching algorithm is ap- 
plied. PSF extraction was completed for both FIRS 1 
and the background field; the final spectrum is the dif- 
ference between them. 

The resulting spectrum from 7.4 to 37.2 /im is shown 
in Figure O Each module is reduced separately, so the 
degree of agreement gives some idea of the reliability of 
calibration. Both the full reduced spectra (solid lines), 
and the average fiux in wavelength bins of AA ~ 0.75 /im, 
1.5 /zm, and 1 /xm for SLl, SH, and LH, respectively, 
are shown. The SH spectrum has the lowest signal to 
noise, so it is binned on the coarsest grid. Only data 
with signal-to-noise greater than one (SLl, SH) or three 
(LH) are included in the binned points. Binned fiuxes 
are listed in Table [TJ 

In Figure [2] the silicate absorption band at 9.7 /im is 
clearly visible, and a hint of the CO2 ice band at 15 /im 
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is also visible. A number of finer features in the LH 
spectrum are most likely real, but will not be discussed 
here. 

3.2. CARMA 230 GHz map 

Continuum observations a,t v — 230 GHz (A = 1.3 mm) 
were completed with the Combined Array for Research 
in Millimeter- wave Astronomy (CARMA), a 15 element 
interferometer consisting of nine 6.1m and six 10.4 m an- 
tennas. Data were obtained in the B-array (100—1000 m 
baselines), C-array (30 — 350 m), D-array (11 — 150 m), 
and E-array (8 — 66 m) configurations between 2007 Oc- 
tober 24 and 2008 December 31. These data were com- 
bined to provide wii-coverage from 4.5 to 500 fcA. Small 
7-pointing mosaics were made in the compact configura- 
tions (D and E) in order to mitigate spatial filtering by 
the interferometer, and to more fully map the spatially 
extended protostellar envelope. 

All three correlator bands were configured for contin- 
uum, 468 MHz bandwidth, observations. A bright quasar 
(1751-1-096) was observed approximately every 15 min- 
utes to be used for complex gain calibration. Absolute 
flux calibration was accomplished using 5 minute obser- 
vations of Uranus, Neptune, or MWC 349. The over- 
all calibration uncertainty is approximately 20%, from 
the reproducibility of the phase calibrator flux on nearby 
days. A passband calibrator, typically 3C454.3, was ob- 
served for 15 minutes during each set of observations, and 
radio pointing was performed every two hours. Observa- 
tions in the most extended B configuration utilized the 
Paired Antenna Calibration System (PACS) to correct 
for phase variations on minute timescales (see L. Perez 
et al. 2009, in preparation). 

Calibration and imaging were accomplished 
with the MIRIAD data reduction package 
(jSault. Teuben. fc Wrightl 119951 ). The resulting 230 GHz 
map of FIRS 1 is shown in Figure [31 with maps made 
at three resolutions: short baseline data only (D and E 
configurations), all data, and long baseline data only (B 
and C configurati ons). The direction of the 3.6 cm jet 
(|Rodriguez et all 119891 : iCuriel et all 11993^ is shown for 
reference. 

The map including all data was inverted with natu- 
ral weighting, cle aned with a Steer CLEAN algorithm 
(jSteer et al.lll984[ ). and restored with an 0.94" x 0.89" 
beam. The rms noise level in the central region is 
6.7 mjy beam^^, the peak and total flux from a Gaus- 
sian fit are 0.42 Jy beam"^ and 1.38 Jy (PA= -4deg), 
and the deconvolved FWHM size is 1.6" x 1.2". The 
synthesized beam corresponds to approximately 240 AU, 
while the longest baselines provide a resolution better 
than 100 AU (0.46" x 0.40"). A gaussian fit to the long 
baseline data yields PA= 25 deg, approximately 75 deg 
from the 3.6 cm jet axis (PA= —50 deg). 

The extended, complex nature of the source is appar- 
ent, thanks to the excellent wii-coverage achieved with 
multiple configurations. Although it is difficult to see the 
more extended envelope even in the short baseline map, 
it is clearly visible as an amplitude peak at itw-distances 
< 20 fcA in a plot of amplitude versus uw-distance (Fig- 
ure S]). Note that the interferometer does filter out flux 
at uf-distances less than 4 fcA, corresponding to the sep- 
aration of the closest antenna pairs. Figure [J shows that 
most of the source flux is concentrated at low and in- 



TABLE 1 

Spitzer IRS AND BROADBAND Spitzer, SHARC II, and 

BOLOCAM FLUXES USED IN THE SED FITS 



Wavelenj 
{^J,m) 


;th Flux 

(Jy) 


Uncertainty 

(Jy) 


Aperture 
(arcsec) 


Instrument 


3.6 


0.00085 


0.00016 


2.2 


IRAC 


4.5 


0.0026 


0.0005 


2.2 


IRAC 


5.8 


0.0023 


0.0005 


2.2 


IRAC 


7.8 


0.004 


0.002 


3.7 


i-tvo 


f'm 1 ^ 


8.5 


0.0018 


0.0010 


3.7 


l-tvo 


?m 1 \ 


9.1 


0.0011 


0.0006 


3.7 


irvo 


?m 1 \ 


9.7 


0.0002 


0.0003 


3.7 


TTJ ^ 

i-tvo 


f'm 1 ^ 

^Oiji ) 


10.3 


0.0016 


0.0008 


3.7 




?m 1 \ 


10.9 


0.0017 


0.0009 


3.7 


irv-O 


?m 1 \ 


11.5 


0.0024 


0.0012 


3.7 


TTJ ^ 

i-tvo 


f'm 1 ^ 


12.1 


0.0034 


0.0017 


3.7 


l-tvo 


?m 1 \ 


12.7 


0.005 


0.003 


3.7 




1 \ 

^Oijl j 


13.3 


0.005 


0.002 


3.7 


i-tvO 




13.9 


0.006 


0.003 


3.7 


IRS 


(SLl) 


14.5 


0.007 


0.003 


3.7 


IRS 


(SLl) 


15.0 


0.0019 


0.0011 


3.7 


IRS 


(SLl) 


15.4 


0.009 


0.006 


10.5 


IRS 


(SH) 


16.3 


0.021 


0.008 


10.5 


IRS 


(SH) 


18.1 


0.033 


0.010 


10.5 


IRS 


(SH) 


19.8 


0.121 


0.017 


11.1 


IRS 


(LH) 


20.7 


0.21 


0.04 


11.1 


IRS 


(LH) 
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IRS 


(LH) 


22.7 


0.59 


0.07 


11.1 


IRS 


(LH) 


23.8 


1.04 


0.11 


11.1 


IRS 


(LH) 


25.1 


2.0 


0.2 


11.1 


IRS 


(LH) 


26.5 


3.1 


0.3 


11.1 


IRS 


(LH) 


28.0 


4.6 


0.5 


11.1 


IRS 


(LH) 


29.7 


6.8 


0.7 


11.1 


IRS 


(LH) 


31.6 


9.1 


1.0 


11.1 


IRS 


(LH) 


33.9 


11.8 


1.3 


11.1 


IRS 


(LH) 


36.0 


13.4 


1.5 


11.1 


IRS 


(LH) 


70 


77 


41 


17 


MIPS 


350 


195 


79 


40 


SHARC II 


1100 


5.9 


1.2 


40 


Bolocam 


Note. - 


— Apertures 


in which fluxes 


are calculated correspond 



to the instrument PSF, with the exception of the SHARC II and 
Bolocam fluxes. Calibrated IRS fluxes are averaged in wavelength 
bins of AA ~ 0.75 /xm, 1 /xm, and 1.5 for SLl, SH, and LH, 
respectively, as described in Section 13.11 The IRS flux is the 
mean within a bin, and the instrumental uncertainty is the error 
in the mean {a/y/N). All uncertainties include a 10% systematic 
uncertainty in addition to the instrumental uncertainty. 

termediate ww-distances (extended structure), but the 
source is clearly detected at ww-distances greater than 
200 fcA, indicating an unresolved or marginally resolved 
compact (< 1") component. Values of the 230 GHz flux 
as a function of uw-distance are given in Table [21 

3.3. Spitzer, Bolocam, and SHARC U broadband data 

Broadband infrared data for FIRS 1 are ta ken from the 
"Cor es to Disks" Spitzer Legacy program ([Evans et al.l 
1200^ . which imaged approximately 1 square degree in 
the c loud with IRAC and MIPS ([Harvev et alJ 120061 . 
|2007| ). The same region was m apped at A = 1.1 rn m with 
the Bolocam bolometer array enn et al.l 120031) at the 
Calte ch Submillimeter Observatory (CSO) ( Enoch et al.l 
|2007[) . These data provide wavelength coverage from 
A = 3.6 to 1100 Aim (IRAC 3.6, 4.5, 5.8, 8.0 /xm; MIPS 
24, 70, 160 /m; Bolocam 1100 ^im). FIRS 1 is not de- 
tected in the 2MASS catalogs. 

Broadband fluxes are used to determine the bolometric 
luminosity and temperature (11.0 Lq and 56 K), and 
are included in the model fits in Section [5l below. The 
total envelope mass (8.0 M©) is calculated from the total 



The Resolved Massive Disk in FIRS 1 



5 




10 -10 -20 ^ „ ^ _ 4 2 -2 -4 

A e ^ 5 0-5 -10 AC ^ 

Arc Seconds . ^ , Arc Seconds 

Arc Seconds 



Center: R.A. 18 29 49.80 Dec +01 15 20.4 

Fig. 3. — CARMA 230 GHz maps of Serpens FIRS 1 for short baseline data only (D,E configurations; panel A), all data (panel B), and 
long baseline data only (B,C configurations; panel C). Contours in panel B are (2, 4. ..10, 15, 20, 30. ..70) times the Icr rms of 6.7 mjy beam~^, 
for a synthesized beam of 0.94" X 0.89" (shown, lower right). Contours in panels A and C are similar but start at Acr and 6a, respectively, and 
pane l A has additi onal contours at (90, 110, 130 cr). Note the change in scale in each panel. The direction of the 3.6 cm jet HHodrieuez et aEI 
119891 ; [Curiel et al...l993) is shown for reference. 
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Fig. 4. — CARMA 230 GHz visibility amplitude versus uii-distance for Serpens FIRS 1. Observations in the B, C, D, and E CARMA 
antenna configurations provide tji'-coverage from approximately 4.5 kX to 500 kX. The expected value in the case of zero signal, or amplitude 
bias, is indicated by a dotted line and is typically small (less than 0.1 Jy). 



flux in a 40" aperture at A = 1.1 mm, assuming the 
envelope is optically thin at 1.1 mm, a dust opacity o f 
nimm = 0.0114 cm g-1 (|Ossenkopf fc Hennindll994D 
and a dust temperature of Tn = 15 K fsee lEnoch et al.l 



l2009l for more details). 

We also include in the observed SED the 350 fim con- 
tinuum flux (M. Dunha m et al. 2009, in prep aration), ob- 
tained with SHARC-II ijDowell et al.|[2003l ) at the CSO. 
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TABLE 2 
CARMA 230 GHz visibilities 



MD-distancc 


Flux 


Unccrtainty 


(kA) 


(Jy) 


(Jy) 


4.50 


2.45 


0.12 


7.50 


2.05 


0.07 


10.5 


1.96 


0.07 


13.5 


1.70 


0.06 


16.5 


1.60 


0.06 


19.5 


1.38 


0.05 


22.5 


1.31 


0.05 


25.5 


1.17 


0.04 


28.5 


1.10 


0.04 


31.5 


1.07 


0.04 


34.5 


0.97 


0.04 


37.5 


0.89 


0.04 


40.5 


0.86 


0.04 


43.5 


0.83 


0.04 


46.5 


0.75 


0.03 


49.5 


0.73 


0.04 


52.5 


0.74 


0.04 


55.5 


0.68 


0.03 


58.5 


0.77 


0.04 


61.5 


0.72 


0.04 


64.5 


0.67 


0.04 


67.5 


0.61 


0.04 


70.5 


0.62 


0.05 


73.5 


0.54 


0.04 


76.5 


0.52 


0.04 


79.5 


0.57 


0.05 


82.5 


0.70 


0.05 


85.5 


0.71 


0.06 


88.5 


0.59 


0.07 


91.5 


0.54 


0.05 


94.5 


0.67 


0.05 


97.5 


0.47 


0.05 


100.5 


0.36 


0.06 


Note. — 


Visibilities and uncer- 



tainties used in the model fits. The 
amplitude bias, or expected value 
for zero signal, has been subtracted 
from the data. [The full table will 
appear in the online version.] 

The SHARC-II flux samples the peak of the SED, and 
helps constrain the long-wavelength side of the model 
SED. All fluxes used in the SED fit are given in Table [H 
including uncertainties, aperture diameters and instru- 
ment used for the observations. 

4. RADIATIVE TRANSFER MODEL 

To model the observed emission from FIRS 1, we use 
the two-dime nsional Monte Carlo rad i ative transfer code 
RADMC of iDullemond fc Domini^ (|200l . RADMC 
performs both Monte Carlo radiative transfer to derive 
the temperature distribution from an input density dis- 
tribution, and ray tracing to produce images and pho- 
tometry in specified apertur es. We adopt a den sity pro- 
file very similar to that of iCrapsi et al.l ()2008[ ) , which 
includes three components: the envelope, the outflow 
cavity, and the disk. For both the envelope and disk 
we use the dust opac i ties f rom Table 1, column 5 of 
lOssenkopf fc Henningl |1994) for dust grains with thin 
ice mantles, including scattering, interpolated onto the 
necessary wavelength grid. We note that although this 
is a young source, there could be significant difference 
in the dust properties of the disk and envelo pe, as have 
been demonstrated in some Class I sources ()Wolf et "all 

[2oba . 



The envelope density profile is that of a rotating, col- 
lapsing sphere ijUlrich et al.lll967f) : 
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where = cos 9, Rc — Rcent is the centrifugal radius, 
and po is the density at {r,6) = (i?c,0), which is set by 
the total envelope mass M^nv = 8 M©. A gas-to-dust 
mass ratio of 100 and mean molecular weight of 2.33 are 
included. Here /ip — cos is the solution of the parabolic 
motion of an infalling particle, given by: 



Mo - M 
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ticent J VMoSm /io. 

The model has no time dependence; Rcent is not a true 
centrifugal radius defined by the conservation of angu- 
lar momentum during collapse, but simply a set radius 
where the density peaks, and inward of which the density 
drops to very low values. The envelope outer radius Rout 
is the maximum value of the radial grid, so the envelope 
density is zero for r > Rout- The outfiow cavity is de- 
fined by setting the density to zero in the region where 
cos 00 > cos(Ang/2). This results in a funnel-shaped 
cavity, which is conical only at large scales where Aug is 
the full opening angle. 

The disk density is given by a power law dependence 
in radius and a Gaussian dependence in height: 



Pdisk{r, 0) 



2nH{r) \Rd 



disk 



exp 



disk 



M, 



disk ■ 



(3) 
At 



where Eg is set by the input 
Rdisk pl changes from —1 to —12, effectively setting 
the disk radius. The scale- height variation (flaring) is 
given by: H{r) = r{Ho/Rd^sk){r/Rd^sky^■ We set 
p2 = 2/7, corresponding t o the self-irradiated passive 
disk of iChiang fc Goldreiclil ()1997D . Given that the disk 
is in large part hidden by the envelope, a much simpler 
description would likely wor k just a s well, but we choose 
to follow the setup of .Crapsi et al.l (12008. ) here. Most of 
the disk parameters are held fixed for the main model 
grid, with only Mdtsk and Rdisk varying. 

Table [3] lists the range of input values for the model 
input parameters, some of which are held fixed. The in- 
ternal luminosity is set by the bolometric luminosity; it 
is input to the model as the stellar luminosity, although 
most likely a majority of the luminosity is due to accre- 
tion. We do not include an interstellar radiation field 
here, as the luminosity of any reasonable field is negligi- 
ble compared to the internal luminosity of FIRS 1, and 
has only a negligible effect on the long- wavelength SED. 

Small single parameter grids were run to test that 
"fixed" parameters have no significant affect on the 
model SED or millimeter visibilities. The scale height 
(Hq) and fiaring (pi) of the disk do affect the 3 — 20 /im 
fiuxes, although only for Incl < 10 deg, where there is too 
much NIR flux to match the data regardless of the value 
of these parameters. A very puffy disk {Hq > 0.5Rdisk) 
produces more emission in this range, while a very thin 
disk {Hq < O.lRdisk) produces less emission. Similarly, 
a flared disk produces more MIR emission than one with 
no flaring. 
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TABLE 3 

Range of parameter, values used in the radiative transfer model grid 



Parameter 


Fixed? 


Range 


Description 




Protostar 


^star 


Y 


11 Lq 


Internal luminosity 


Tstar 


Y 


4000 K 


Protostar effective temperature 






Enve 


lope and outflow 


Menv 


Y 


8.0 Mq 


Total mass of envelope 


Rout 


N 


3000 - 12000 AU 


Outer radius of envelope 


Rcent 


N 


50 - 1000 AU 


Centrifugal (inner) radius of envelope 


Ang 


N 


5-80 deg 


Outflow opening angle 


Incl 


N 


5-90 deg 


Inclination angle 


Disk 


Mdiak 


N 


0.0-3.0 Mq 


Disk mass 


Rdisk 


N 


50 - 1000 AU 


Disk radius 


Ho 


Y 


0.2 B.disk 


Disk vertical pressure scale height 


pi 


Y 


-1.0 


Disk surface density radial power law (r < Rdisk) 


p2 


Y 


2/7 


Power law for H(R) (disk flaring) 



Note. — The internal luminosity is set by the bolometric luminosity of the source, determined 
fro m the broadba nd SED, and the envelope mass is set by the 1.1 mm Bolocam single dish flux 
fsee lEnoch et al7^2009). Ang is the full outflow opening angle. Incl is the line of sight inclination 
angle of the disk: dog is face-on, 90 deg is edge-on. Stellar, envelope, and disk parameters are 
discussed in Section |4] 



5. RESULTS 

To determine the best fit envelope parameters we run 
a grid of models varying Rout, Rcent, Ang, and Incl. 
This results in a total of 588 envelope models, each 
observed at 13 inclination angles. A nominal disk of 
Mdisk = 0.01 Mq, Rdisk = 150 AU is used for aU en- 
velope models. A separate grid varying Mdisk and Rdisk 
with fixed envelope parameters includes 140 models. Our 
tests show that the disk has little effect on the SED for 
this source, making separate grids feasible. 

The model grids are compared to the SED and 230 
GHz visibilities with a analysis using data from Ta- 
bles [1] and m Contour plots of the resulting reduced 
(x^) are shown in Figure [5l where contours for both the 
SED and 230 GHz visibilities are plotted. The model 
and data visibilities are calculated by the same method, 
using vector averaging in radial annuli. 

The best-fitting model (Rout = 5000 AU, Rcent = 
600 AU, Ang = 20 deg, Incl = 15 deg, Md^sk = 1-0 Mg, 
Rdisk — 300 AU) is compared to the data in Figure [6l 
The same envelope model with no disk is also shown for 
reference, as is an envelope model with no inner envelope 
hole {Rcent = 20 AU; dotted line). Determination of the 
best fit model from Figure [5] is described in Sections 15.11 
and 15.21 Due to the small uncertainties and the limited 
number of models, the reduced values for even the 
best-fit model are still fairly high (x^ ~ 13 and 5 for the 
SED and visibilities, respectively). The bolometric tem- 
perature and luminosity of the best-fit model are 35.3 K 
and 17.3 Lq. The bolometric luminosity differs from the 
input stellar luminosity due to inclination effects. 

While we only show the best-fit model here, there is a 
range of values for each parameter that can reasonably fit 
the data, as determined by eye from x^ plots and visual 
inspection of the SED fits. We find that the envelope 
parameters have the following reasonable ranges: Rout ^ 
5000-7000 AU, Rcent ~ 400-600 AU, Ang ~ 10-30 deg, 
and Incl ~ 10 — 25 deg. Reasonable disk parameters are 
Md^sk ~ 0.7 - 1.5 Mq, and Rd.sk ~ 200 - 500 AU. 

Literature fluxes are shown in comparison to our data 
and the best-fit model SED in Figure [71 Shown are IRAS 



HIRES 25, 60, and 100 fim from lHurt fc Barsonvl (I1996H 
SCUB A 450 and 850 /im peak fiuxes from Davis et al.l 
(IT999I). JCMT 800 1100, 1300, and 2000 fluxes from 
Casah etal] (119931) . a nd the OVRO 3 mm flux from 
Testi &: SargentI ()1998[ ). They cannot be compared di- 



rectly to the model SED because many of them are peak 
fluxes calculated in small apertures; circles show the 
corresponding model values computed in similar aper- 
tures. While there is not perfect agreement, the model 
is roughly consistent with the literature values, with the 
exception of IRAS fluxes. The IRAS observations are 
lower resolution than the Spitzer maps, and FIRS 1 may 
be confused with nearby protostars (24 — 45" away). Lit- 
erature fluxes are not included in the x^ fitting. 

5.1. Envelope Structure 

Envelope parameters are determined first, using the 
grid of envelope models (Figure [5] (A)-(C)). The 4-D x^ 
space is collapsed along the parameters not plotted in 
each panel. We conclude below that Rout < 7000 AU 
based on the 230 GHz visibility x^ contours in panel 
(A), so only models with Rout < 7000 AU are included 
in panels (B) and (C) to avoid complicating the plots 
(because many models with Rout > 7000 AU can fit the 
SED). 

With a few exceptions, the SED is much more sensi- 
tive to envelope parameters than are the 230 GHz visi- 
bilities. Rout, however, is only mildly constrained by the 
SED and is somewhat degenerate with Rcent- Increasing 
either parameter lowers the total opacity of the envelope, 
allowing more MIR emission to escape. In addition, for 
Rout > 7000 AU, the SED provides no constraint on 
Rcent because the opacity through the envelope is al- 
ready quite low. In this case, the 230 GHz data do help 
constrain the envelope parameters because visibility am- 
plitudes at uw-distance < 30 kX trace extended emission. 
A narrow peak in the uv plane (Figure |4|) corresponds 
to a large envelope outer radius, and a wider peak to a 
small outer radius. The best compromise between the 
visibilities preferring smaller Rout and the SED prefer- 
ring larger Rout is Rout 5000 AU (Figure [5] A) . The 
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Fig. 5. — Reduced contours resulting from fitting the observed SED and CARMA 230 GHz visibilities (Tables [T] and [2]l to the grid 
of envelope and disk models. Models have been run for the full parameter ranges shown; the model grid resolution corresponds to the 
axis labels (e.g. Rcent values of 50, 100, 200, 300. ..1000 AU), but the distribution has been smoothed for a better visual representation. 
Contours from fits to the both the SED (magenta) and visibilities (cyan, tick marks indicate downhill direction) are shown, although 
envelope parameters (panels A-C) are primarily constrained by the SED, while disk parameters (panel D) are constrained by the millimeter 
visibilities. In panels (A)-(C) the distribution is collapsed along the parameters not plotted. The lowest contour plotted is x^ = 8. 



value of Rcent witli the lowest SED in this case is 
600 AU. 

A more compact envelope creates a high opacity to 
shorter wavelength emission, thus a larger Rcent is 
needed to decrease opacity close to the protostar and 
to match the observed NIR and MIR emission. For 
Rout < 7000 AU, the range of reasonable centrifugal radii 
are Rcent = 400 — 700 AU. There are a few models with 
small Rcent and low x^; when the viewing angle is just 
down the edge of the outflow (Incl = Ang/2) in a dense 
envelope {Rout — 4000 AU). In these special cases the 
opacity is lowered just enough to give a similar emergent 
SED as models with a large inner envelope hole. 

A compact envelope with Rout ~ 5000 AU is consis- 
tent with the crowded star formation region in which 
FIRS 1 is located. The nearest embedded protostar that 
is known to have an envelope is 45", or approximately 
11000 AU away, and seve ral other embedded sources are 
within a few arcminutes (|Enoch et al.ll200"9[ ). While we 



do not know the actual 3-D distances, envelopes with 
radii Rout ~ 5000 AU are certainly reasonable in this 
clustered region. 

Both the outflow opening angle and inclination are 
fairly narrowly constrained by the SED (Figure [5] 
(B),(C)). There is a degeneracy between Aug and Incl; 
all models where the line of sight is directly within the 
outflow cavity (Incl < Ang/2) have very high val- 
ues because they produce a large excess of NIR emis- 
sion. The SED is best flt by models with low inclinations 
(Incl < 35 deg) ; larger inclinations produce very high ex- 
tinction in the MIR and cannot match the observed MIR 
flux. Low inclinations may be in conflict with the orien- 
tation of the 3.6 cm jet, which has been interpreted as 
being almost in the plane of the sky based on proper mo- 
tion of emission knots in the jet (Moscadelli et al. 2006). 

To summarize, the short uu-spacing visibility data fa- 
vor small outer radii, while the SED favors larger outer 
radii, with the best compromise lying at Rout ~ 5000 AU. 
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Fig. 6. — Best-fit envelope and disk model, compared to the observed SED (left) and 230 GHz visibilities (right). The SED includes the 
binned IRS spectra and broadband data from 3.6 — 1100 fj,m (Table[T]l. Three models are shown for comparison: the best-fit envelope and 
disk (red), the best fit envelope with no disk (blue), and a model with no inner envelope hole {Rcent = 20 AU; dotted). The envelope-only 
model {Mdisk = Mq) indicates the relative contributions of the envelope and disk to the 230 GHz visibility amplitudes. 



For Rout = 5000 AU, the SED is minimized for 
i?^g„j = 600 AU (panel (A) of Figure O . With Rcent set 
to 600 AU, it is straightforward to determine the best-fit 
Ang and Incl from the SED curves in panels (B) and 
(C). 

None of the models accurately reproduce the shape of 
the silicate absorption feature or the slope of the spec- 
trum from 10 — 20/im, although the uncertainties on the 
observed fluxes are also quite high in this region. This 
is most likely a feature of the dust model and not the 
envelope density profile. Similarly, the dust model does 
not include the 15 /xm CO2 ice absorption feature. 

Note that the input stellar luminosity (11 L©) is 
based on the bolometric luminosity calculated before 
the SHARC II 350 fim was available, and thus is lower 
than the bolometric luminosity of the best-fitting model 
(18 L0). Increasing the stellar luminosity to 18 Lq does 
not change the results dramatically; the best-fit centrifu- 
gal radius is a bit lower, 400 AU, without changing the 
other parameters. The reasonable range of Rcent is much 
larger however, allowing for Rcent as low as 50 AU for 
larger outflow opening angles (e.g. 40deg). 

5.2. Disk Structure 

After the envelope parameters {Rcent, Rout, Ang, Incl) 
have been determined, we run a separate grid in disk 
mass and radius with envelope parameters fixed. The 
resulting contours are shown in Figure [5] (D). Disk 
parameters are entirely constrained by the CARMA 
230 GHz visibilities; the SED is insensitive to both Mdisk 
and Rdisk- A quite massive {Mdisk ^ 0.7 — 1.5 M©) 
and resolved {Rdisk ^ 200 — 500 AU, compared to the 
maximum resolution of 100 AU) disk is required to ac- 



count for the significant flux at intermediate uw-distances 
(20 - 100 fcA; Figure 0. 

Typically, only a lower limit can be placed on the disk 
mass because once the disk emission becomes optically 
thick larger masses do not increase the millimeter flux. 
Here, however, because the disk is resolved the mass is 
more tightly constrained. For FIRS 1, fitting the visi- 
bilities out to maximum ww-distances from 50 — 300 kA 
produces the same best-fit disk parameters, although 
the goodness of fit increases with more data. The uv- 
coverage required to get a good fit should depend on the 
disk and envelope structure. 

Disk properties derived by fitting the 230 GHz visibil- 
ities are relatively insensitive to the assumed envelope 
parameters. It is true that a somewhat less massive 
disk would be required if there were no inner envelope 
hole; for example for Rcent = 50 AU, the best-fit is for 
Mdisk ~ 0.4, Rdisk ~ 400 AU. However, for envelope pa- 
rameters with any reasonable fit to the SED (for exam- 
ple, Rcent = 400 AU, Rout = 7000 AU), Md^sk = 1.0 M© 
and Rdisk — 300 AU remain the best fit to the 230 GHz 
visibilities. 

The disk mass is also reasonably robust to uncertain- 
ties in the envelope mass (8 M©) and 230 GHz cali- 
bration. For Menv = 10 M© the best-fit disk is un- 
changed, because only the amplitude of the narrow peak 
at small uv-distances is affected by the envelope mass. 
For Menv = 6 Mq, a slightly larger, more massive disk 
{Mdisk ^1-5 M©, Rdisk ~ 400 AU) is required to account 
for the decrease in envelope fiux at small uv distances. 
The overall fit is poorer than for Menv = 8 M©, however. 
Systematic uncertainties in the CARMA 230 GHz fluxes 
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Fig. 7. — Best-fit model compared with literatu re fluxes: IRAS 
HIRES 25, 60, 100 /im fHurt & BarsomJ [ToM) . SCUBA 450, 
850 Atm ( peak fluxes; D avis et al. 1999), JCMT 800, 1100, 1300, 
2000 fim (Casali et al. 1993), and OVRO 3 mm IITesti fc Sargent! 
[1998). The model has been computed in apertures consistent with 
the literature measurements (circles). IRAS fluxes are based on 
lower resolution data than the Spitzer fluxes used here, and may 
be confused with nearby protostars. Literature fluxes are not used 
in the SED fit. 

have a shghtly larger effect, with a 30% change in overaU 
flux calibration producing a corresponding 30% change in 
the disk mass: Mdisk ^ 0.5 — 0.7 M© for a 30% decrease, 
and Mdisk ^ 1-5 — 2.0 Mq for a 30% increase. 

5.3. Other models 

Here we compare our models to other disk and enve- 
lope models. This serves both as a check on our derived 
envelope and disk parameters by indicating which con- 
clusions are dependent on the density model, and a test 
t hat our models give the m ost reasonable fit to the data. 

iHo gerheiide et al.l (|1999l ) observed FIRS 1 with the 
OVRO interferometer at 3.4, 2.7, and 1.4 mm. They used 
a power law envelope model plus a point source, with 
the dust temperature power law set at —0.4 and fixed in- 
ner and outer envelope radii of 100 and 8000 AU. These 
millimeter data were best fit by an envelope with mass 
6 Mq and density power law —2.0, plus an unresolved 
point source of approximately 0.7 Mq. IHogerheiide et ahl 
()1999f l note that they are unable to separate the inner en- 
velope from any disk emission, and thus cannot place a 
meaningful limit o n the disk itself. Compared to the 
IHogerheiide et all (|1999[ ) ■uii-coverage, 10 — 180 kX at 
1.4 mm, our CARMA observations trace much more of 
the envelope (down to 4.5 fcA), allowing us to separately 
model and remove the envelope contribution. Our data 
also have much higher signal-to-noise on long baselines 
(uii-distances > 100 fcA). 

In addition to the rotating, collapsing spheroid (or "Ul- 
rich" ) envelope models described in Section lU we also 



ran a small grid of models with a simple power law en- 
velope density profile (p cx r^^), plus a conical cavity. A 
steep power law, p cx r~^, provides a reasonable fit to 
the visibilities without requiring a massive disk or large 
inner envelope hole. The best fit is for Mdisk = 0.1 M©, 
Rdisk = 100 AU, and Rcent = 100 AU, as shown in Fig- 
ure [51 Here the emission at intermediate ww-distances 
is filled in by the envelope, which reaches high densi- 
ties close to the protostar, thus requiring less disk emis- 
sion. But, only a special combination of outflow opening 
angle and inclination can come close to matching the 
SED (Ang = 40deg, Incl = 25deg or Ang = 60deg, 
Incl — 30deg; looking down the edge of a wide outflow 
cavity) , and even the best-fit model gives a much poorer 
flt to the SED than the Ulrich models. In general, the 
power law models seem unable to reasonably fit the NIR 
and MIR emission, although only p ~ —2 and —1.5 have 
been tested here. 

The ability of the power law envelope model to fit the 
visibilities without a large disk is consistent with some 
previous studies that have found that disks are often not 
req uired to mill imeter data of Class and Class I sources 
(e.g iLoonev et a l. 2000). The results here demonstrate 
that it is necessary to include both spectral and visibility 
data in order to fit a consistent disk and envelope model. 

We use the online SED fitting tool of IRobitaille et al] 
(|2007f ) as another estimate of the envelope parameters, 
although the 230 GHz vis ibiliti es cannot be included in 
the fit. [Robitaille et all (|2007f ) use an envelope setup 
similar to ours, with the envelope infall rate Me„„ set- 
ting the fiducial density po (rather than Mgnv as used 
here). The best- fit model corresponds to a protostar 
with age t = 2 X 10^ yr, = 1.8 Mq, = 7 R©, 

= 4400 K, envelope infall rate Menv = 10"^ M© yr~\ 
= 11000 AU, Ang = 27deg, and Incl = 75deg. The 
total luminosity and envelope mass of this best flt model 
are consistent with our values (18.4 L© and 7.4 M©). 

Looking at the 10 best fltting models, only the age 
(i < 2 X 10^ yr), protostellar mass (M* < 2 M©), tem- 
perature (r* = 3000 — 4500 K), and envelope infall rate 
{Menv — 10"^ — 10""* M© yr~^) are reasonably well con- 
strained, while the other parameters cover a large range. 
We do not attempt to constrain the disk properties as 
the SED is relatively insensitive to the disk in embedded 
sources. In addition, we cannot use the online grid to 
constrain the envelope inner radius, because it is fixed 
to the disk inner radius and very few models with both 
large envelope mass (> 1 M©) and large inner envelope 
radius (> 10 AU) are included in the grid. This differ- 
ence in the inner envelope behavior likely accounts for 
the large outer radius and inclination required by the 
IRobitaille eFaLl (|2007[ ) models. 

Given our limited exploration of various models, we 
feel that the Ulrich envelope model provides the best 
description of the observed SED. While the derived disk 
parameters do depend on the input envelope density pro- 
file, even in the most conservative case Mdisk ^0.1 M©. 

6. DISCUSSION 

Our derived disk mass of Mdisk ~ 1.0 M© within a ra. - 
dius of 300 AU is consistent with the [Brown et alj (l2000f) 
limit of Mdisk > 0.1 M©, as well as the lHogerheiide et al.l 
()199a) limit of 0.7 M© on the unresolved mass within 
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Fig. 8. — Best-fit power law envelope model, p oc r~^, compared to the observed SED and 230 GHz visibilities. A steep power law 
alleviates the need for a massive disk {M^nsk =0.1 Mq, Rdisk = 100 AU is the best fit), but the power law models are unable to match 
the observed MIR Spitzer fiuxes. 



a radius of 100 AU. The early evolutionary state of 
FIRS 1 is confirmed by the low bolometric temperature 
{Tboi ^ 56 K) and the small disk-to-envelope mass ratio 
{Mdisk/Menv ~ 0.1) despite the high disk mass. Thus, 
our results suggest that large disks can accumulate very 
early in the protostellar collapse process. 

The FIRS 1 disk is likely too small to be considered 
a magnetically supported "pseudo-disk" . Given the ex- 
pected young age of FIRS 1, however, both the mass 
and radius derived here are much larger than expected 
for di sk formation via grav itational collapse of a rotating 
core. iTerebev et al.l ()1984[ ) predict that the disk radius, 
where centrifugal balance is achieved, should depend on 
the initial rotation f2 and isothermal sound speed Cg in 
the core as: 

" ^ (o.35km s'O (4XIG-14S-1) (lov) 

(4) 

Based on the statistical relationsh ip between Tboi and 
time derived in lEnoch et"all (|2009[ ) {Ttoi oc t^-^], the 
bolometric temperature of FIRS 1 suggests that it has 
an age of 0.7 — 0.8 x 10^ yr. For a reasonable sound 
speed, Cs ~ 0.23 km s~^ (T ^ 15 K), this age and 
Rdisk — 300 AU requires an initial rotation rate of ap- 
proximately 5 X 10"^^ s~^. This value is higher than 
typical dense cores, wh ich have fl ^ 10^^^ — lO"^** s~^ 
(| Goodman et al.iri993f) . Alternatively, if the age of the 
source is actually closer to 3 x 10^ yr, a more typical 
rotation rate would be sufficient for growing a 300 AU 
disk. 

Similarly, we can estimate how long it would take for 
the disk to build up 1 Mq via infall from the envelope. 
For an infaU rate of Menv ^ c^/G ^ 10~^ Mq yr~\ and 



conservatively assuming that all of the infalling material 
falls onto the disk rather than directly onto the proto- 
star, accumulating 1 Mq would take 10^ yr. Although 
this is probably close to the age of FIRS 1, a disk of 

1 Mq requires that little of the infalling material be ac- 
creted from the disk onto the star. Below we mention 
a few plausible methods for building up a large circum- 
protostellar disk in this object. 

A Class lifetime longer than a few times 10^ yr, i.e., 
longer than the estimated timescale for Cl ass sources 
in nearby low mass st ar forming regions (|Enoch et al.l 
\200^ lEvans et al.ll2009[ ) , would allow larger disks to grow 
before the end of the Class phase. 

FIRS 1 may have a higher envelope infall rate than 
average, allowing the disk to quickly accumulate mass. 
The bolometric luminosity of FIRS 1 is quite large 
compared to the typical luminosity of YSOs in nearby 
low-mass star forming regions (< ILq; iDunham et al.l 
,200_8; Eno ch et al. 2009). A luminosity of 11 Lq im- 
plies an accretion rate onto the protostar of at least 

2 X 10-5 Mq yr-i (for ~ R^Ltoi/GM^, i?, - 5 i?© 
and Af, ~ 1 Mq). If this corresponds to an even higher 
envelope infall rate, Mdisk ~ 1 Mq could easily be 
achieved in 0.7 — 0.8 x 10^ yr. 

If the FIRS 1 disk has very low viscosity, mass may 
build up in the disk with very little accreting onto the 
protosta r (although this is at odds with the high lumi- 
nosity). iBrown et al.l (|200C1| ) suggest that early disk for- 
mation and similar disk masses in the Class and Class I 
phases could be achieved with a time-dependent viscos- 
ity, low at early times and higher by Class I. 

Perhaps more likely, this source could have recently en- 
tered a period of relatively rapid ac cretion, as expected in 
the episodic accretion scenario (e.g lHartmann fc KenvonI 
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[1985L lEnoch et all [200l lEvans et all [2001 . Such a 
high mask disk around a presumably low protostar mass 
should be unstable, and undergoing rapid accretion, ex- 
plaining the high luminosity. In this picture, the current 
high accretion phase would have been preceded by a pe- 
riod of low accretion onto the protostar while material 
built up in the disk (assuming infall from the envelope 
onto the disk is steady). 

A larger sample is certainly needed to determine if 
such large, high mass disks are typical in the Class 
phase. The recent J0rgcnsen ct al. (2009) study of 20 
Class and Class I sources finds disks masses in Class 
from 0.01 - 0.5 Mq, and Mdrsk/Menv ratios of 1 - 10%. 
If we calculate o ur di sk mass by the same method as 
iJorgensen et al.l ()2009D . which uses the flux at 50 kA 
and assumes an optically thin, unresolved disk, we get 
Mdisk = 0.6 Mp. Wh ile this is at the high end of the 
iJorgensen et al.l (|2009| ) disk sample, the disk to envelope 
mass ratio (~ 8% when using M^isk = 0.6 M0) is con- 
sistent with their results, as is the idea that disks are 
already well established in the Class stage. 

Regarding the envelope structure, it is important to 
keep in mind that while we refer to Rcent as the cen- 
trifugal radius, our model is not dynamical, and there 
is no dependence on rotation rate. Thus the sharp drop 
in density inside of this radius could have any number of 
causes, including a companion that has cleared out mate- 
rial, as well as rotational collapse onto a disk. Any binary 
companion with a disk mass larger than 0.1 Mq should 
have been detected. There is a tentative second detec- 
tion 500 AU to the northwest of FIRS 1 (see Figure[3]C); 
the peak of approximately 55 mJy would correspond to a 
disk mass of ~ 0.06 M©, but this may just be a "clumpy" 
feature in the envelope. The more likely explanation is 
that the inner envelope cavity is a result of collapse in 
a rotating core, and creation of the 300 AU radius disk. 
Alternatively, a clumpy envelope could cause a similar 
decrease in MIR opacity and might alleviate the need for 
an inner envelope hole ijlndebetouw et aLll2006f) . 

If the disk and envelope are physically connected, with 
both the disk and inner envelope hole governed by ro- 
tation in the collapsing core, we might expect Rdisk ~ 
Rcent- Although the best-fit Rcent is a factor of two larger 
than Rdisk here, the range of reasonable values allow for 
Rcent to be as small as 400 AU and Rdisk to be as large 
as 500 AU. Thus, a physical continuity between the disk 
and envelope is certainly plausible. 

7. SUMMARY 

We utilize Spitzer IRS spectra, high resolution 
CARMA 230 GHz continuum data, and broadband pho- 
tometry together with a grid of radiative transfer models 
to characterize the disk and envelope structure of the 
Class protostar Serpens FIRS 1. Our conclusions are: 

1. Radiative transfer models combined with mid- 
infrared spectra and millimeter data with excellent uv- 
coverage can reasonably constrain envelope parameters, 
including the inner (centrifugal) radius, outer radius, and 
outflow opening angle. In all cases there is a range of 
parameter values able to reasonably fit the data. Once 
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the envelope parameters have been determined, the mass 
and radius of the disk are robustly constrained by mil- 
limeter interferometry data with ww-coverage from < 5 
to > 300 kX. 

2. We find a centrifugal radius for FIRS 1 in the range 
400 — 700 AU, indicating a large "hole" in the inner en - 
velope, similar to IRAS 16293 (jJorgensen et al.ll2005[ ). 
Unlike IRAS 16293, however, there is no strong evidence 
for a binary companion that might have cleared out the 
inner envelope. Other explanations for such a large Rcent 
in this source are: (a) collapse of the inner envelope onto 
the disk due to the conservation of angular momentum 
in a rotating, collapsing core, or (b) the Rcent does not 
indicate a true inner radius, but rather a "clumpy" enve- 
lope with much lower opacity in the MIR than a smooth 
envelope density profile. 

3. Using envelope parameters set by the SED, the 
CARMA 230 GHz visibilities require a quite massive, 
resolved disk. The best-fitting disk has a mass of IMq, 
and a radius of approximately 3 00 AU. While this mass is 
consistent with pr evious limits (jHogerheijde et al1ll999l : 
iBrown et al.ll2000h . it also indicates that protostars can 
accumulate relatively massive disks at very early times. 
This is somewhat at odds with theoretic al expectations 
that disks start small and grow with time (jTerebev et al.l 
[1981 . The range of reasonable disk and envelope param- 
eters does allow for a physical continuity between the disk 
and inner envelope. 

4. Our results for FIRS 1 demonstrate the feasibility of 
using this method to characterize the disk and envelope 
structure in a larger sample of Class sources. Similar 
modeling for an unbiased sample will allow us to char- 
acterize the typical disk mass, size, and inner envelope 
structure during the Class phase. 
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